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Abstract. The central notion of this work is that of a functor between categories of 
finitely presented modules over so-called computable rings, i.e. rings R where one can 
algorithmically solve inhomogeneous linear equations with coefficients in R. The paper 
describes a way allowing one to realize such functors, e.g. HoniR, <S)r, Ext^, Torf , as a 
mathematical object in a computer algebra system. Once this is achieved, one can compose 
and derive functors and even iterate this process without the need of any specific knowledge 
of these functors. These ideas are realized in the ring independent package homalg. It 
is designed to extend any computer algebra software implementing the arithmetics of a 
computable ring R, as soon as the latter contains algorithms to solve inhomogeneous linear 
equations with coefficients in R. Beside explaining how this suffices, the paper describes 
the nature of the extensions provided by homalg. 



1. Introduction 

In the setup of finitely presented module categories, homalg realizes functors as mathe- 
matical objects which, up to now, can be composed and derived. To this end it realizes, 
unlike all present systems, the functors not only by the object part, but automatically also 
by the morphism part. 

homalg is abstractly designed and therefore can be used as an extension of other math- 
ematical software providing the necessary ring arithmetics in any concrete problem. 

1.1. homalg as a programming environment for homological algebra. Functors 
map objects of a source category to objects of a target category and, in a compatible way, 
morphisms between two objects in the source category to morphisms between their images 
in the target category. So when one implements a functor one has not only to take care of 
how it acts on objects, but also of how it acts on morphisms between these objects. 

Homological algebraic constructions [CE99, HS97, MR01, Rot 79, Wei94] are present in 
most of the computer algebra systems such as Macaulay 2 [GS] , Singular [GPS05]/Plural 
[LS03], and CoCbA [CoC]. One can often find a procedure, let us call it Horn, to com- 
pute Hom R (A, B) for two finitely presented modules A and B over commutative rings 
that are implemented in these systems. Now given two further modules M and N and 
a morphism M N, what about applying Hom R (A, — ) or Hohir(— to the mor- 
phism ipl Mathematically it is clear how Hohir(A, — ) or Hohir(— , B) induces a morphism 

Rom R (A,M) UomR{Aip \ R omR (A,N) or Uom R (N, B) UomR{ip ' B \ Hom R (M,E), but on 



Date: 2007. 



2 



MOHAMED BARAKAT & DANIEL ROBERTZ 



the level of a computer implementation the information to compute these induced mor- 
phisms is normally not contained in the procedure Horn. So one needs to write or find 
two completely independent procedures, say HomMap and Hom2Map implementing the com- 
putations of the induced morphisms. HomMap is for example essential 1 if one wants to 
compute Ext R (A, B) := R l Rom R (—, B)(A), the i th right derived functor of Eom R (-,B) 
applied to A. Again, a procedure, let us call it Ext, implementing for all % the computa- 
tion of the extension modules Ext^(A, B) is normally not of much help when it comes to 
computing the induced morphisms Ext^(<£>, B) or Exh R (A,(p). Writing the corresponding 
procedures, say ExtMap and Ext2Map, requires one to know about the derivation process, 
and of course to have Horn, HomMap and Hom2Map predefined. So what about composing 
functors? Consider Ext J R (Ext^(-, -), -) for example (cf. [Roo62], [Bj679, p. 58ff]). For 
the object part one needs to compose the procedure Ext with itself using the first argu- 
ment. Let us call this compositum ExtExt. The morphism part has now three components. 
Let us consider Ext^(Ext^(A, —),B) for example. Again, a procedure ExtExt is offhand 
not of much help. For this morphism part one rather needs a predefined Ext2Map to first 
compute ip := ~Ext R (A,tp), then a predefined ExtMap to compute Ext 3 R (ip, B). Let us call 
the resulting procedure ExtExt2Map. 

So in order to derive or compose functors, one needs their part on objects but also their 
part on morphisms. This means, that for each functor one has to implement one procedure 
for the object part and as many procedures as needed for the morphism part. For more 
complex functors, constructed out of given ones via iterated compositions and derivations 2 
this quickly becomes unfeasable. So on the level of the computer implementation the 
following became unavoidable: 

• Include the mathematical information of how the bifunctor Hohir(— , — ) acts on 
morphisms inside the procedure Horn itself. This has to be done in a way, that a 
general procedure, say FunctorMap, is able to extract this information out of Horn. 
Further HomMap and Hom2Map should now be defined only using FunctorMap applied 
to Horn. 

• Implement a general right derivation procedure for a contravariant functor (given 
by both its parts), let us call it RightDerivedCof unctor. Now define Ext only 
using RightDerivedCof unctor applied to the Hom-functor given by both 3 its parts 
[Horn, [HomMap, Hom2Map]]. The procedure Ext, or rather the derivation procedure 
used to define Ext, should be able to reconstruct the mathematical information of 
how the bifunctor Ext R (—, — ) acts on morphisms, and this in such a way that the 
same procedure FunctorMap mentioned above is able to extract this information out 
of Ext. The derivation procedure should reconstruct this information alone from 



Maybe not in full generality, since one only needs to apply Hom^(— ,B) to morphisms between free 
modules R lxli ^R lxl *. 

2 Left derivation for covariant functors and right derivation for contravariant functors (=cofunctors). 

3 Since we right derive HomR(— ,— ) with respect to its first argument, the corresponding morphism 
procedure HomMap has to be marked in the input of RightDerivedCof unctor. This is a minor technical 
issue. 
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its input [Horn, [HomMap, Hom2Map]]. ExtMap and Ext2Map should now be denned only 
using FunctorMap applied to Ext. 
• Implement a general composition procedure for two functors (given by both their 
parts), let us call it ComposeFunctors. Then define ExtExt only using Compose- 
Functors applied 4 to [Ext, [ExtMap, Ext 2Map]] and [Ext, [ExtMap, Ext2Map]]. The 
procedure ExtExt, or rather the composition procedure used to define ExtExt, 
should be able to reconstruct the mathematical information of how the trifunctor 
Ext^.(Ext^.(— , — ), — ) acts on morphisms, in such a way that the same procedure 
FunctorMap as above is able to extract this information out of ExtExt. The compo- 
sition procedure should reconstruct this information, alone out of its input, which 
is here [Ext, [ExtMap, Ext 2Map]] and [Ext, [ExtMap, Ext 2Map]]. Again, ExtExtMap, 
ExtExt2Map and ExtExt3Map should now be defined only using FunctorMap ap- 
plied to ExtExt. 

Like FunctorMap, the procedures RightDerivedCof unctor and ComposeFunctors should 
be implemented in a way that is independent of the functors they are applied to. So, with 
the same RightDerivedCof unctor one should be able to define K Ext j R (Ext k R (A, -) , B) 
and again with the same FunctorMap compute its part on morphisms, etc. 

Starting from Section 5 we will try to isolate the mathematical ideas that helped us 
to realize functors as mathematical objects. That defining more complex functors in both 
their parts is now an easy, even automatic task is the first defining property of homalg. 

1.2. homalg as a meta-package. Not a single algorithm to compute in any sort of rings 
is implemented in homalg. Rather it is a package built up of homological algebraic defi- 
nitions and constructions. On the highest levels one finds the construction of connecting 
homomorphisms and long exact sequences, the processes of composition and derivation of 
functors and definitions of various specific functors. On lower levels homalg goes all the 
way down till it reaches two procedures, which are basically the only ones required from 
any software implementing the ring specific arithmetics of the (not necessarily commuta- 
tive) ring R. In what follows we refer by u ring package" to such software. To describe the 
two procedures let K be a (left) i?-submodule of the free module R lxn given by finitely 
many generators. Further let b be an arbitrary element of R lxn and A an i?-submodule 
of K, again given by finitely many generators (A < K < R lxn 3 b). Let M £ /£ fcxn resp. 
N £ R axn be the matrix having as rows the generators of K resp. of A. 

(Z) The procedure DecideZero(b, M): 

Effectively decide if b is an element of K or not. "Effectively" means: In case 
the element b belongs to K, then the procedure DecideZero returns zero and, if 
asked to, is able to express b as an i?-linear combination of the generators of M. 
Otherwise some element b' £ b + K is returned. This is equivalent to deciding the 
solvability of the inhomogeneous i?-linear system xM = b and, if asked to, finding a 



Since we compose Ext#(— , — ) (here with itself) with respect to the first argument, the corresponding 
morphism ExtMap has to be marked in the first argument of ComposeFunctors. 
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particular solution x G R lxk , in case one exists. This is simply the straightforward 
generalization of the ideal membership problem to submodules. 
(S) The procedure SyzygiesGenerators(M, N): 

Compute a generating set of the -R-module of solutions of the homogeneous i?-linear 
system xM = mod N. For z G R lxn the statement z = mod N means that there 
exists a y G R lxa , such that z = yN. One calls every such solution a syzygy among 
the generators of K modulo A. 
As the reader may have noticed, it is not required in (Z) that the output b' of DecideZero 
only depends on b + K, i.e. homalg does not require from DecideZero to provide a normal 
form modulo K for every element b G R lxn , except for b G K, where DecideZero must 
return zero. Deciding if two elements b, c G R lxn represent the same class modulo K, is 
reduced to checking if b — c G K, i.e. if DecideZero (b — c, M) = 0. 

In practice however one often solves (Z) by first constructing a different set of generators 
for K out of the given one, where this new set satisfies certain properties 5 , such that a 
reduction algorithm with respect to a set of generators having these properties is available. 
In what follows we call such a set of generators a basis (this should not be confused with a 
free basis). Normally, such a basis also provides a way to algorithmically solve (S). homalg 
enables the ring package to specify a procedure called BasisOf Module to compute such a 
basis. Internally, homalg only uses the procedure DecideZero to perform reductions with 
respect to such a basis. In case the ring package performs reductions without computing 
such a basis the procedure BasisOf Module is to be set to the identity procedure. In what 
follows we will refer to the output b' of DecideZero(b, M) as the reduction of b modulo K 
(or modulo M). 

Summing up one can say: 

(1) homalg is designed to be easily extendable by sparing the user the technical details 
of homological constructions. 

(2) homalg is a meta-package that is designed to easily extend mathematical soft- 
ware implementing ring arithmetics. All what homalg needs from such an im- 
plementation are two procedures: One to effectively solve the ideal membership 
problem (Z) and one to compute a generating set of syzygies, i.e. to solve (S) 
(cf. [Sch80, CL092, KR00, BGTV03, DL06]). We will call such rings computable. 

We stress the following point: For homalg it is irrelevant how (Z) and (S) are solved. 
Hence, the irrelevance of explaining, or even mentioning how to solve these two problems for 
every particular class of rings is a defining property of homalg, namely its second defining 
property. 

In Section 2 we describe the categories homalg is dealing with. In Section 3 we deduce 
from the two procedures DecideZero and SyzygiesGenerators all the other procedures 
used in the sequel. After introducing some notation in Section 4 we describe homalg's 
philosophy of implementing functors in Sections 5 and 6, whereas Section 7 describes how 
homalg derives functors. Appendix A outlines the various ring packages that have been 



5 E.g. a GAUSSian triangular basis in case of fields or an involutivc or Grobner basis in case of poly- 
nomial rings. 
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successfully used in connection with homalg. Appendix B includes some comments to the 
current implementation. Finally, examples are given in Appendix C. All appendices are 
available at the site of homalg [BR07c]. 

It goes without saying that we will suppress some technical issues of the package for the 
sake of mathematical clarity. The best (technical) guide to the package remains its source 
code, which in nearly all parts can be read as a mathematical text. 



2. Presentations 

Let R be a left noetherian ring with one. Denote by C mod the category of finitely 
generated left -R-modules, which is a full subcategory of the abelian category of left R- 
modules. And denote by C mat the category of finite /^-presentations with objects being 
finite dimensional matrices over R, where one identifies two matrices M G ^ lX '° and M' G 
firixio w itu the same number / of columns to one object, if R lxll Vi = i? lxri M', as R- 
submodules of R lxl °. The set Morgmat (M, L) of morphisms between two objects M G R llXl ° 
and L G R l 'i xl 'o is the set of all / X 1' i?-matrices ip with R lxll M(p < R lxl '^L, where one 
identifies two matrices ipi and tf2 to one morphism, if they induce the same i?-module 
homomorphism from coker(i? lxil —> R lxl °) = R lxl ° / R lxll l / [ to coker(i? lx ''i — > R lxl 'o) = 
R lxl 'o/R lxl 'iL. 

As usual a presentation is given by generators and relations. If one takes the classes 
of the Iq standard basis vectors of R lxl ° as generators of M := coker(.R 1>dl R lxl °) = 
R lxl ° j R lxll M, then the h rows of M are the defining relations between these generators. We 
thus call M a presentation matrix or a relation matrix for M. Further we call R lxll ¥l = im(M) 
the relation subspace for M generated by (the rows of) M and denote it by (M). 

If M contains a unit in the j th column, then the row containing this unit is a relation 
expressing the j th generator of M — coker(M) as an i?-linear combination in terms of the 
other generators. Hence one can, using elementary matrix transformations, rewrite the 
matrix of relations with respect to the remaining generators and obtain a new relation 
matrix with one less column (and one less row). One iterates this process until the relation 
matrix is free of units. So, without loss of generality, one can assume that the relation 
matrix M is free of units. For deciding invertability of a ring element and computing its 
inverse homalg uses the procedure Leftinverse described in 3.1.3. 

Summing up, the functor coker that maps 

R hxh 3HmM:= coker( J R lx/l ^ R lxl °) 

is an equivalence between the categories C mat and C mod . Its inverse functor is to fix a 
presentation matrix M G Objg ma t for each module M G 0bj emO d and to use the generators 
of the fixed presentations to express maps in Mor e mod as matrices in Morgmat (cf. [GP02, 
p. 101]). 

In what follows we no longer distinguish between the two categories and therefore we 
simply denote both by C. 
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3. Basic procedures 
3.1. Procedures based on DecideZero. 

3.1.1. RightDivide(B, A, L). This procedure is the ring-theoretic version of finding a partic- 
ular solution of an inhomogeneous linear system of equations. Hence, it is not astonishing 
that in a lot of computations it plays the decisive role (cf. 3.1.2-3.1.4). 

For three i?-matrices A, B and L with the same number of columns, one wants to compute 
an X with 

B = XA mod L, 
i.e. an X, such that there exists a Y with 



B = XA + YL = (X, Y) 




After computing (X, Y) one often simply throws away Y. So without loss of generality we 
forget about L, i.e. consider the system B = XA. To find X one computes a basis G for A 
(i.e. a matrix G with rows being a basis for the rows of A) together with a matrix Ca, such 
that G = CaA. Then one reduces B modulo G, i.e. computes a matrix N := DecideZero (B, G) 
(cf. 1.2,(Z)), with rows being those of B reduced modulo G. We compute N together with 
a matrix C B , such that N = B + C B G. If N is not the zero matrix, then the system is not 
solvable. Otherwise X = — C B C A is a solution. Thinking of X as "X = BA _1 " justifies the 
name of the procedure. 

The most prominent application of RightDivide in homalg is the following situation: 
Given the three modules and the two morphisms 



M' 




N 



satisfying the so-called image condition im7 < im/5 (as submodules of N), find a third 
morphism ip : M' — > N' that makes the triangle commute. The image condition is obviously 
a necessary condition for such a ip to exist. There are two instances of special importance 
for us in which such a ip always exists: 

(1) M' is a free module (of finite rank r), or 

(2) (3 is injective. 

In the first case let (pi, . . . ,b r ) be a free basis of M' . Define biip to be any element of 
the set (6j7)/3 _1 , i.e. any element in the preimage of b^ under (3, and since M' is free on 
(6i,..., b n ) this extends by linearity to a morphism ip : M' — > N' satisfying %p(3 = 7. In the 
second case (3 is an isomorphism onto im/3 and one defines ip := 7/3 -1 , where 7 is viewed 
morphism M' — > im /3. 

Now let M', N' and N be relation matrices of M', N' and ./V respectively. Then we are looking 
for a matrix ip, together with matrices Y and Z such that 

(i) 7 = ^/3 + ™, 
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(ii) M> = ZW. 

If one regards ip merely as a matrix, i.e. not necessarily a morphism M' — > N', then 
condition (i) is nothing else but the above mentioned image condition. Therefore, for the 
situation described above, such a matrix 7 always exists. But now condition (ii) states 
that if; carries the relations of M' to relations of N', i.e. it requires ip to be a morphism 
from M' to N' . While condition (i) is already in the form needed by RightDivide 6 , ip 
on the left hand side of condition (ii) is multiplied from the right 7 . It now happens that 
in the two instances (1) and (2) condition (ii) is automatically fulfilled, whenever (i) is. 
In the first instance one needs to additionally require that M' is given on a free basis. 
Then the matrix of relations M' is the zero matrix (with one row and r columns), so (ii) is 
trivially satisfied for any ip by taking Z to be the zero matrix of appropriate dimension. 
Case (2) is equivalent to saying that the preimage of the subspace (N) of relations for N 
(generated by N) under j3 coincides with the subspace (N') of relations for N' (generated 
by N'): ((N))/?" 1 = (N'). But since 7 is a morphism M' — > N any relation m! G (M') of M 

satisfies = m'7 = (m'ip)j3 mod (N). Hence, there exists an n' e (N'), such that n' = m'tp. 
So ip takes relations of M' to relations of N' and is hence a morphism M' — > N' . 

3.1.2. CompletelmSq. This is the most prominent reincarnation of RightDivide in homalg. 
We call an incomplete square of the form 

•p 

P 

N N 

an image square if the image condition im(aip) < im(/3) (as submodules of N) is satisfied. 
This is precisely the above situation for 7 = oap. In Section 5 CompletelmSq is applied to 
image squares with injective j3 (case (2)) and in Subsection 7.2 to ones where M' is free 
(case (1)). In these two instances, as shown above, the image square is completable by a 
morphism ip : M' N' which is directly computable using RightDivide. 

3.1.3. Lef tinverse. The typical applications of Leftinverse correspond to the two situ- 
ations (1) and (2) described for RightDivide in 3.1.1 (cf. Example C.3): 

(1) Either take M' = N free (case 3.1.1,(1)), 7 = id and (3 : N' — > N surjective. Then 
the image condition is trivially satisfied. 1/) : M' = N — > N' is then nothing else 
but the left inverse or a split of j3. 



Following the notation used in defining RightDivide we set A = /?, B = 7 and L = N. 
7 If the ring R is commutative, one can use the Kronecker product to construct for any matrix C over 
R a matrix C satisfying row(CA") = row(X)C for any composablc matrix X, where row(C) is the row 
vector consisting of the rows of C written behind each other in the obvious order. By this trick one can 
rewrite the two conditions (i) and (ii) in a single affine condition, where row( , 0), row(y) and row(Z) are 
multiplied from the left as required by RightDivide. Since the resulting affine system is, in general, much 
bigger compared with the initial ones solving it is computationally expensive. Cf. [ZL02]. 
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(2) Or take M' = N, 7 = id and j3 : N' — > N an isomorphism. Then in particular /5 
is injective (case 3.1.1,(2)), and the image condition is trivially satisfied, ip : M' = 
N — > N' is then nothing else but the (left)inverse of f3. 

a 

3.1.4. Preimage. Let M — > L be a morphism between the two modules M and L := 
coker(L). If the matrix B consists of rows that are in the image of A, i.e. if the image 
condition is satisfied, then the rows of the matrix X computed by RightDivide are the 
preimages of the rows of B under A. In Subsection 7.4 we use the procedure Preimage to 
construct connecting homomorphisms. 

3.2. Procedures additionally based on SyzygiesGenerators. For a matrix A £ R llXl ° 
a matrix of generating syzygies or a syzygy matrix is a matrix X £ R l ? xl ^ _R lx/2 X = 

ker(i? lx ' x — > R lxl{> ). A slight generalization of this is when one requires that R lxl ' 2 X = 
ker(R lxl1 coker(L)), where L £ R axl ° has the same number l of columns as A. One then 
calls X a syzygy matrix of A modulo L, or relative to the module coker(L). 

So the rows of the matrix X of generating syzygies generate the solution space of the 
homogeneous linear system xA = mod L. For z £ R lxl ° the statement z = mod L 
means that there exists a y £ R lxa , such that z = yL. 

3.2.1. ResolutionOfModule(M). By iterating the process of taking syzygies one obtains a 
free resolution of the module M = coker(M): 

. . . Lpl+2 } fi}xli+i filxk ... ^ Rlxh yi=M ) _Rlx/ m _^ 

of desired length. What homalg additionally does is the following: As explained in Section 2 
we can assume the relation matrix (or the first syzygy matrix) M = </?i free of units. Starting 
from i — 2, whenever the i th syzygy matrix ipi is computed, homalg uses the units appearing 
in it to locate the redundant rows of the (i — l) st syzygy matrix This means, if <pi 

contains a unit in the j th column, then the row in ifi containing this unit says that the j th 
row of ipi-i is an i?-linear combination of the other rows (of ifii-i), so it is redundant for 
generating the (i — l) st syzygies and can be omitted from </?i-i- The j th row is thus deleted 
from (fii-i and the 2 th syzygy matrix tpi is recomputed. The rows of the latter are again 
used to locate further redundant rows of </?i-i- This process obviously stabilizes. Then one 
proceeds with the (i + l) st syzygy matrix ifi+i, etc. We end up having a free resolution 
where all the matrices </?j are free of units. 

In the case of graded modules over a positively graded ring with degree part a field, 
or modules over local rings, this process indeed yields a minimal free resolution (cf. [Eis95, 
p. 472] and [Sch03, p. 40]). In both cases eliminating units from the matrices tfi of the 
resolution is equivalent to requiring that all entries belong to the unique maximal graded 
resp. unique maximal ideal of R. 

However, and although in general this process does not yield a "minimal" 8 free resolution, 
it often enough reduces the involved dimensions of the matrices considerably (cf. [KR05]). 



Such a notion is in general not even well defined, cf. [Eis95, p. 472]. 
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3.2.2. Subf actorModule(Mi, M 2 ). Given two matrices Mi and M 2 with the same number of 
columns, computing a presentation matrix for the subfactor module ((Mi) + (M 2 ))/(M 2 ) 
(short-hand: (Mi)/ (M 2 )) goes as follows: One computes a basis B of M 2 . Then one reduces 
Mi modulo B and gets N. The syzygy matrix S of M modulo N is the desired presentation 
matrix. 



4. Categories of complexes of given finite length 

As always, let C be the module category of Section 2. By D k (G) we denote the category 
of chain complexes C : Ck — > Ck-i C\ of C of finite length h — 1 and their chain 

maps. 

Definition 4.1 (The category D\{G\). For i £ {1, . . . , k} let D k (G) be the factor category 
of D k (G) defined by forgetting in a chain map the morphisms between all but the i th chain 
module. 



Note that for a morphism in D\ (C) 



a 



a 



i+l 



i-l 



a 



a 



i-l 



there exists at least one commutative completion 



CI 



k 



i+l ~ 



Ch 



a 



i+l 



a 



c. 



-Ci 



c. 



i-l 



which is a preimage in D h (Q). 

By definition D\(G) is simply a different notation for C. From now on all functors we 
consider take values in the module category G. 



Definition 4.2 (The functor Objf ) 
a complex to its i th module: 



Denote by Gbjf : D k (G) -> G the full functor mapping 



a 



Ck 



i+l 



a 



i+l 



a 



C' 



i-l 



a 



i-l 



C"; 



Definition 4.3 (The functor Morf ). 
a complex to its (i + l) st morphism: 



Denote by Morf : Df (G) 



C\ Ci 
— ► C the full functor mapping 



CI 



Ck 



^i+l 



n ^+1 n 
^i+1 *" 



a 



i-l 



Ci 



W+l *" 
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5. The morphism part of a functor 

A functor is by definition a map between two categories that maps objects to objects 
and, in a compatible way, morphisms to morphisms. To be allowed to speak about functors 
one needs to be able, not only to define the functor on objects, but also on morphisms 
between objects. 

It might be puzzling for the reader that we start with a section devoted to the morphism 
part rather than the object part of functors. The reason for that will become clear towards 
the end of this section. 

Here we write everything for covariant functors. Adapting things for contravariant func- 
tors is done in the obvious way. 

For a morphism S — » T and a functor F we want to compute F(<p). There are three 
cases homalg distinguishes: 

(Cmp) The functor F is defined in homalg as a composition of two functors: F = F\o F2. 
(Der) The functor F is defined in homalg as the i th (left) derivation 9 of another functor: 
F = UC 

(Bsc) The functor F is defined in homalg neither by composition nor by derivation. In 
what follows such functors are called homalg- basic 10 functors. 

Dealing with (Cmp), i.e. with a composed functor is easy: Since F(ip) = Fi(F2(<p)) one 
is able to reduce computing F(<p) to computing -F 2 (y2) and then Fi(F 2 ((^)). A bit more 
involved is the case (Der), i.e. when F = Lj G. In Subsection 7.2 it is shown how to reduce 
the computation of F(ip) to essentially computing G(tp). For homalg-basic functors the 
idea is to reduce the computation of F{ip) to completing an image square (cf. 3.1.2, case 
(2)). To this end one embeds F(S) resp. F(T) in a module HuhV(S') resp. HuhV(T) with 
an induced morphism Hu11f(<£>) determined by <p: 

(Hull) F(S) — HuhV(S') 



F{<p) 



-F(T)— Hull F (T) 

Expressed categorically, HuhV is a functor, which we call a hull functor of F, and lf is a 
natural transformation, which we call the corresponding natural embedding. 

For a given functor F the idea is to either define the hull functor HuhV from scratch (e.g. 
for 11 Cokemel (in 6.1.1), HomJR. (in 6.2.1), TensorProduct (in 6.2.2) and Horn (in 6.2.3)) or 



homalg provides procedures to left derive covariant functors and right derive contravariant functors. 
These are the cases computed via a projective resolution of the module. 

10 The prefix homalg indicates that the notion of "basic functor" is not a mathematical definition. If, 
for example, the functor Homij(HorriR(— , R), R) would have been implemented without using homalg's 
composition procedure ComposeFunctors (see B.4), we then would call it homalg-basic. 

^Since one can directly provide the action of Cokernel and TensorProduct on morphisms, one takes 
them as their own hull functors. 
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to define Hull^ by composition of already defined functors and the forgetful functors Obj^ 
(Def. 4.2) and MorJ- (Def. 4.3) (e.g. for Kernel (in 6.1.2) and DefectOf Horns (in 6.1.3)). 

The procedure in homalg that precisely accomplishes the above mentioned distinction is 
called FunctorMap (cf. Appendix B.7). 

6. Some homalg-BASic functors 

Since we need to explain how to construct for each functor its hull functor, we need to 
first of all mention which specific standard method for computing the object part we use 
(see for example [GP02]) and then how to construct the hull functor in the specific setup. 
The hull functor is the missing piece of data that allows one to automatically compute the 
functor part on morphisms. 

All functors we consider are of the form F : Df(C) — > C. We simply distinguish two 
types of functors, depending on whether F{C) is a subfactor module of Obj^(C), for all 
objects C G -Df(C), or not. In what follows a functor of the previous kind is called a 
functor producing subfactor modules. 

To avoid confusing the two parts of a functor with source category 12 -D 2 (C) (which has 
as its set of objects also morphisms between modules), we use two different names for the 
two parts: If the object part of a functor is called F, then its morphism part is called FMap. 
The name F is also used to refer to the functor in both its parts. 

6.1. Functors producing subfactor modules. 

6.1.1. The functor Cokernel. On the level of objects the covariant cokernel functor asso- 
ciates to a morphism between two modules its cokernel. Cokernel : D\(Q) — > C, i.e. 

A' - > A Cokernel (a) 



CokernelMap((/j) 



B'^-^B Cokernel ((3) 

such that 

A 1 — A »- Cokernel (a) ^ 

<fi CokernelMap(</3) 

B' *- B *■ Cokernel (/3) 

is commutative and exact. 

Defining the object part Cokernel (a) is simple: After fixing a presentation matrix A for 
A and generators for A' one can view a as a matrix with the same number of columns as 

A. Now take the supermatrix ( ^ J as a presentation matrix for Cokernel (a). 

The cokernel functor is the most basic functor in our setting, in the sense that com- 
puting its morphism part CokernelMap(^) is trivial: If we take the residue classes of the 
generators of A resp. B to be the generators of Cokernel (a) resp. Cokernel (/?), then 



12 



This is the case in 6.1.1, 6.1.2. 
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the matrix representing CokernelMap((/?) coincides with that of (p. Hence one can set the 
hull functor Hull Co kernei = Cokernel and the natural embedding tcokemei is the identity 
transformation 13 . 



6.1.2. The functor Kernel. On the level of objects the covariant kernel functor associates 
to a morphism between two modules its kernel. Kernel : -Df(C) — > C, i.e. 



.4 



B 



A" 



B" 



Kernel(ct) 

| KernelMap((p) 

Kernel^) 



such that 







Kernel(a) 



KernelMap^) 

^ Kernel^) 



A 



B 



A" 



B" 



is commutative and exact. 

Defining the object part Kernel (a) goes like this: After fixing a presentation matrix 
A resp. A" for A resp. A" one can view a as a matrix with the same number of columns 
as A". Compute the syzygy matrix iota of a modulo A" (cf. the beginning of Subsection 
3.2). Now use Subf actorModule to compute Kernel(a) := (iota)/(A) (cf. 3.2.2). See also 
[GP02, p. 97]. 

In the case of the functor Kernel the hull functor is obviously Hull Ke rnei = Obji? an d the 
natural embedding t Kernel is the transformation embedding the kernel of a morphism into 
its source module: 



0- 



-Kernelfa 



^Kernel (<*) 



■A = 0bj2(A A A") 



The matrix of the natural embedding is simply iota. 

6.1.3. The functor DefectOf Horns. On the level of objects the covariant defect functor 
associates to two composable morphisms ct\ and with a^ai = their defect of exactness 
ker(a 2 )/im(ai). DefectOf Horns : D\{£) — > C, i.e. 



A' 



B' 



A 



A" 



B 



B" 



DefectOf Horns («!, a 2 ) 

DefectOf HomsMap(</3) 

Def e ct Of Horns (fii, fa) 



This is not the whole truth. The functor Cokernel calls a procedure named Presentation, which, 
among other things, tries with the help of the procedure BetterGenerators, to reduce the number of 
generators of the resulting module using either normal form algorithms referred to in the introduction and 
at the end of Appendix A, in case they are available and applicable, or otherwise, beside the one described 
in Section 2, several clever heuristics. Taking care of this, is a simple, but technical issue. 
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such that 



*- Def ectOf Homs(o:i, a 2 ) *■ Cokernel(ai) " 2 > A" 



Def ectOfHomsMap(<£>) 



CokernelMap (ip) 



02 

»- Def ectOf Horns (ft, ft) > Cokernel(ft) B" 

is commutative and exact. 

The definition of the object part of Def ectOfHoms uses the ideas in 6.1.1 and 6.1.2: Fix 
a presentation matrix A for A and A" for A". Fix generators for A'. Then one can view 
cki resp. a 2 as a matrix with the same number of columns as A resp. A". Using A, a 2 and 
A" compute the matrix iota of the embedding of Kernel(a 2 ) in A as in 6.1.2. Then use 

Subf actorModule to compute the defect Def ectOf Homs(ai, a 2 ) ■= (iota)/ /( a * 

In the case of the functor Def ectOfHoms the hull functor is obviously Huli Defect0fH oms = 
Cokernel o Mor^ and the natural embedding ^DefectofHoms is the transformation embedding 
the defect of two composable morphisms into the cokernel of the first morphism: 

„ ™ .... / s 'DefectDf Horns (ai,a2) , s 

s-Def ectOfHoms (ai, a 2 ) ^Cokernel (cti), 

where Cokernel(ct!) = Cokernel (Mor| (A' — > A ^> A")). The matrix of the natural em- 
bedding is iota. 

The functors Cokernel and Kernel are obviously special cases of this functor. 



6.2. Other types of homalg-basic functors. Contrary to functors producing subfactor 
modules, where the morphism part of the hull functor is given by the induced morphism 
on the factor module, here we need to explicitly mention how the morphism part of the 
hull functor is defined. 

For the rest of the subsection let 

M = coker(M) = R lxl ° /R lxh M, L = coker(L) = R lxl '» /R lxl '^L, 

be two finitely presented modules: 

i? 1 ^ 1 — i? lxfe — M -0, Rixii-^Rixi'o^L -0. 

Further let M — > N be a morphism, between two finitely presented modules M and 

N = coker(N) = R lxl o/R lxl h, 

i.e. 

<p G R loXl °. 

Note that we write the morphisms on the right, i.e. we use the row convention. 
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6.2.1. The functor Hom_R. Let R be a (not necessarily commutative) ring with a fixed 
involution 14 9. For a right .R-module H, define the left .R-module H e by setting H e = H 
as abelian groups and r ■ h := h9(r) for all r G R and h G H . So the involution 9 allows 
one to rewrite any right .R-module H as a left R- module H e . 

Recall that Hohir(— , R) is a contravariant functor from the category of left -R-modules 
to the category of right -R-modules. We use the involution 9 to transform the resulting 
right -R-module into a left .R-module. We call the resulting functor Hom^(- , R). The idea 
is to reduce the computation of the homomorphism module to computing a kernel. Since 
Hom fl (- ,R) is left exact, one obtains the exact sequence of right -R-modules: 

Hohir(M, R)^^ Rom R (R lxl °,R)^+Rom R (R lxl \R). 

We compute Hohir(M, R) as the kernel of the right most morphism. Following the row 
convention we identify B.om R (R 1X: ' , R) with B? xl , which justifies the notation (M-) for the 
right most morphism. Applying 9 yields the exact sequence of left -R-modules: 

-Hom^(M, R)^R lxl °^R lxh , 

where (K 9 ) ab = 9(M ba ) and (R lx1 ) 9 = R lxl . Thus set 

Hom_R(M) := Kernel(-R lxio ^ R lxh ). 

This definition proposes setting Hull Ho ni_R(-^0 := -R lx '° and taking the embedding of the 
kernel Hom_R(M) in R lxl ° as the natural embedding (cf. 6.1.2). The morphism part of the 
hull functor is hence defined by Hull Hom _R(v2) := f e : 

^ Hom_R(M) ^ R lxl ° 

HomMap_R((/p) 

> Hom_R(iV) ^ R lxl 'o 

where cp e : R lxl ° — > R lxl 'o is a morphism of free left modules. 

Although the definition of Hull H0 m.R(^) depends on the presentation of M, it is never- 
theless functorial because of the equivalence of the categories C mod and C mat . 

Hom_R is a contravariant functor. The transposition in the definition of M e (and ip e ) is 
the manifestation of this contravariance. 

6.2.2. The functor TensorProduct. Let .R be a commutative 15 ring. Recall that — ®r— is a 
bifunctor, covariant in both arguments. The idea is to reduce the computation of the tensor 
product module to computing a cokernel. Since — ®r — is right exact in both arguments, 
the tensor product of the two presentations R lxl1 R lxl ° of M and i? lx ''i — > R lxl 'o of L 



= (R loXl Y 
V 

= (R i ° xi y, 



An order 2 anti-automorphism: 9 2 = idfj and 9(ab) = 9(b) 0(a) for all a, b e R. 
'Cf. Subsection 6.2.4. 
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(each regarded as a two term complex) is a presentation of M® R L (cf. [GP02, Section 2.7]): 



(R 



lxZi 



R lxl 'o 



1R 



(R 



lxl 



R lxl'^ 



© 



R lxl 



R lxl' 



M 



1R 



0. 



We compute M ® R L as the cokernel of r. After identifying R lxj ® R R lxk with R lx i k the 

M <g> I,/ 
it, ® L 



morphism r is given by the matrix T 



, where ® is the Kronecker product 



of matrices and i/ resp. Iy is the identity matrix on M resp. L. 

Our convention in defining the Kronecker product of two matrices A = (a^) and B 
is the usual one: A ® B := (a^B). 

We define 

TensorProduct (M,L) := Cokernel^ 1 ^ 1 ^ ^ ^ # lxW «). 

Define the morphism part of the functor TensorProduct (—, L) to be the KRONECKER 
product 

TensorProduct (ip, L) := tp (8) Iy . 

Here the hull functor Hull TensorProduct (_ ) L) coincides with the functor and the natural em- 
bedding is the identity transformation. 

Analogously for the functor TensorProduct (L, — ) define 

TensorProduct (L, tp) := ij/ Cg> tp. 

Again, the hull functor Hull TensorProduct (_ i / ) coincides with the functor and the natural 
embedding is the identity transformation. 

TensorProduct (— , — ) is a bifunctor, covariant in each argument. 



is a 



6.2.3. The functor Horn. Let R be a commutative ring. Recall that Hohir(— , 
bifunctor, contravariant in its first argument and covariant in the second. The idea is 
to reduce the computation of the homomorphism module to computing a kernel. Since 
Hom^(— ,L) is left exact for any module L and Hohir(P, — ) is exact for P projective (or 
free) one obtains (cf. [GP02, p. 104]): 





t 







Hom H (M,L) 







Hom iJ (M, R lxl °) 



Hohir(M, R 



lxZi 



^Hom i? ( J R lx '°,L) - 
Rom R {R lxl ° , R lxl '°) 

-\ 

Rom R (R lxl ° , R lxl 'i) 



Hom R ( J R lx ' 1 ,L) 



I'L 



■H.om R (R lxl \R lxl ° 
Rom R (R lxh ,R lxl \ 



We compute Hom R (M, L) as the kernel of k in the first row. Identifying Hom^i? 1 ^, R lxk ) 
with R> xh justifies the notation used for the morphisms of the lower right square of the 
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Cf. Subsection 6.2.4. 
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above diagram. Further, identifying R? xk with R lx i k (by writing all the j rows as one 
long row) gives rise to the identification 17 of HoniR(i? lxi °, L) with Cokernel(J; ® L) and 
of Homfl(i? lxil , L) with Cokernel(/; 1 <g> L). The induced morphism 

K : Eom R (R lxl '\ L) -> Hom fl (i2 lxil , L) 

is then given by the matrix (M ® h' ) tT ■ Thus define 

(M®/,/ ) tr 

Hom(M, L) : = Kernel (Cokernel(/ Zo ® L) 2— > Cokernel(J il ® L)). 

This definition proposes setting Hull Hom( _ jL )(M) := Cokernel(/ i() <g> L) and taking the 
embedding of the kernel Hom(M, L) in Cokernel(/; ®L) as the natural embedding (cf. 6.1.2). 
The morphism part of the hull functor is hence defined by Hull HO ni(-,L)(v 9 ) := (v 9 ® ^i' ) tr '- 

>■ Hom(M, L) Cokernel(/ /o ® L) 

HomMap(v5) '' "(V®^) tr 



»- Hom(iV, L) *- Cokernel(J^ <g> L). 

To address the functoriality of Hom(— , — ) in the second argument, interchange the role 
of M and L, and set Hull Ho m(L,-)(M) := Cokernel(/^ <g> M) and Hull H om(L,-)(v 9 ) := h' ® V 9: 

*- Hom(M, L) »- Cokernel(/ / / ) ® M) 



Hom2Map(t/?) 



> Hom(iV, L) ^ Cokernel(/ i / ® N). 

Hom(— , — ) is a bifunctor, contravariant in its first and covariant in its second argument. 
Transposing the matrix M £g> If (and tp (g> fy) is the manifestation of the contravariance in 
the first argument. 

6.2.4. One last word on the commutativity of the ring R. For the most general definition 
of tensor product of modules one starts with a not necessarily commutative ring R, a 
right i?-module M R , and a left .R-module rN. If the i?-module structure of M resp. 
N comes from a (Q, _R)-bimodule resp. an (R, 5')-bimodule structure, then their tensor 
product over R is in a natural way a (Q, 5')-bimodule qMr <S)r rNs- Q and S, again, are 
not necessarily commutative rings. Analogously, let M and N be two left modules over 
a not necessarily commutative ring R and denote by Hohir(#M, r N) the abelian group 
of i?-homomorphisms. If the i?-module structure of M resp. N comes from an (R, Q)- 
bimodule resp. an (R, S^-bimodule structure, then Hom^(^M, rN) is again in a natural 
way a (Q, 5')-bimodule. 

Note that in both cases the resulting module might not be finitely generated as a (Q, S)- 
bimodule, even if M and N are finitely generated as -R-modules. 



_^ R jxi' becomes jjixjli Ii^ RWo and R l ° xk ^ R hxk becomes R lxl ° k {mIkr ; 
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In the special case where either M or N is an (R, i?)-bimodule, then M N resp. 
Hohir(M, N) is again an i?-module 18 . This is always the case when the ring is commutative, 
since then every i?-module is an (R, i?)-bimodule in the obvious way. 

The Maple implementation of homalg does neither support changing the ring, nor bi- 
module structures. These are issues we want to address in future implementations. But 
even though, the above mentioned problem of non-finite generation will remain the major 
obstacle. 

7. Derived functors 

The philosophy of derived categories is, roughly speaking, to replace a module by one of 
its resolutions, and then to look at the resolution as a very special type of complexes, with 
homology concentrated at degree 0. After inverting quasi-isomorphisms, one obtains the 
derived category, where the objects are quasi-isomorphism types of complexes. Especially, 
all resolutions of a module become isomorphic objects in the derived category. 

Using the cylinder-cone-translation construction [GM03, III. 3] one constructs out of 
every short exact sequence of complexes, a so called distinguished triangle. By passing 
to homology we again obtain distinguished triangles in the category of graded objects 
(which one can view as cyclic complexes, i.e. complexes with zero boundary maps [GM03, 
III.2.3]). A popular way to start, is to take a distinguished triangle coming from a short 
exact sequence of complexes that are simultaneously resolving a short exact sequence of 
modules. Then one applies a functor, that turns such distinguished triangles again into 
distinguished triangles, and at last one takes the homology. The classical way of writing 
such a distinguished triangle of homologies is as a long exact homology sequence. 

The reason for recalling the standard definitions in the following subsections is not only 
to indicate how they are computed using homalg, but to finish the discussion of Section 5. 
This is done in Subsection 7.2. 

7.1. The procedure ResolveModule: Resolve a module. By resolving a module M, 
which we view as a complex concentrated in degree 0, we obtain a complex of free (resp. 
projective) modules and a quasi-isomorphism 



P2- 






\ 




\ 


■o — 


-^o — 


- M 



After inverting quasi-isomorphisms all resolutions become isomorphic. 
For an additive functor F and a module M with a resolution 

P '■ Pq+l *~Pq *~Pq-l " • ■ *~Pl ^-Po ^0 

define the q th derived functor L q F applied to M by setting L q F(M) := H q (F(P)), the 
defect of the two consecutive morphisms F(P q+ i) — > F(P q ) and F(P q ) — > F(P q -i) (cf. 6.1.3 
DefectOf Horns). 

18 In theory, this special case could work in homalg, and indeed it does for very special cases, of course, 
beside the trivial case when M = R or N = R. 
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There is a cheaper method to compute the left derivation of a right exact covariant 
functor, which is based on [HS97, the definition of L q in IV.(lO.l), p. 156 and Prop. IV. 5. 5, 
p. 133] and uses Kernel instead of DefectOf Horns. By duality, there is a cheaper method to 
compute the right derivation of a left exact contravariant functor, which is based on [HS97, 
Prop. IV.5.8] and uses Cokernel instead of DefectOf Horns. Both methods are implemented 
in homalg. 

7.2. The procedure ResolutionOf Seq: Resolve a morphism. Given a morphism 
M N one resolves M and N freely 

(Lift) 



Pq+l 




-P q -i^-- 


■^Pl- 


-Po- 


- M — 


^0 
















P' 


->. P' 

q 




■^P[~ 


P' 


-N — 


-o, 



and computes the tfq's by iteratively completing image squares. Here one needs a free 

resolution of M to be in case (2) of CompletelmSq, 3.1.2. Applying F one gets F(P q ) F ^ Pq \ 
F(P'). But now the object part of the functor L q F applied to M (resp. N) is a subfactor 
module of F(P q ) (resp. F(P')) and one is in the situation of DefectOf Horns, 6.1.3. This 
is all what FunctorMap needs to compute the morphism part of a derived functor. This 
finishes the discussion of Section 5. 



nmsnes tne aiscussion or section o. 

7.3. The procedure ResolveShortExactSeq: Resolving a short exact sequence 

modules. To resolve a short exact sequence of modules — > M' — > M — > M" — > 0, < 



starts with a resolution of M": 



of 

one 




\ 

M' 
\ 

M ■ 
\ 



P" 


- P" - 


_^ pi' 


M" — 


^0 


\ 


\ 


\ 


\ 



















Then one completes the middle line by taking free hulls of iterated pullbacks. Finally one 
fills the upper line by taking kernels to obtain: 

(M) 

\ \ \ \ 

















\ 


\ 


\ 


1 


- p' 


+P[- 


p' 


M' — * 


\ 


\ 


\ 


\ 


-^-P 2 - 


^Px- 






\ 


\ 


\ 


1 


- P" 


- P" ~ 


pll 


-*■ M" — »- 


\ 


\ 


\ 


\ 















0. 
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with exact columns and rows. This method is implemented in the procedure ResolveShort- 
ExactSeq. There are several other methods to resolve a short exact sequence simultane- 
ously, cf. [HS97, Proof of Theorem IV. 6.1]. 

Computing the pullback A' of B' — > B <¥- A is reduced to computing a kernel, namely 

-p) 

that of A © B' > B. The two maps B' <— A 1 — > A of the pullback are the two parts 

(a ip ) 

of the kernel embedding -> A' A® B' (cf. 6.1.2). 

7.4. The procedure LongExactHomologySeq: Connecting homomorphism and long 
exact sequences. Applying an additive covariant functor F to the truncation of the 
diagram (M) 

(P) 














{ 


\ 


{ 




- P' 


-P{- 


pi 


^0 


\ 




\ 






-Pl- 


-P - 


^0 


\ 


* 


\ 




- P" 


- p? - 


P" 


^0 


\ 


\ 


\ 
















results in a diagram P = F(P) where the columns are still exact, but where the rows are 
now in positive degrees, in general, no longer exact. Connecting the homologies of the rows 
one obtains the long exact homology sequence 

■■■^L q F(M') -> L q F(M) - L q F(M") % h q _ x F(M') - • • • , 



where the connecting homomorphisms are computed via the snake lemma applied to the 
diagrams 










I 


P' IB" 

r q+ll n q+l 




\_ 




Pq+l/ Bq+1 


-z q 


\ 


\ 


P" IB" 


- 7" 


\ 










by taking kernels and cokernels, where Bi := im(P i+1 — > Pi) and Zj := ker(P, — > 
Pj-i). One computes the connecting homomorphism from ker(P (?+1 /P (J+1 — > Z q ) to 
coker(P +1 /P +1 — > Z ) by a diagram chase, which accounts in taking preimages twice, 
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namely under the maps P q+ \/B q+ i — ► P" q+l / E>" q+1 and Z q — > Z g . For this the procedure 
Pre image from 3.1.4 is used. 

Applying a contravariant functor to (P) yields a long exact cohomology sequence. The 
corresponding procedure is called LongExactCohomologySeq. 

Appendix A. The ring packages 

The following Maple ring packages have successfully been used with homalg. In each 
of the following descriptions we append a list of rings which can be dealt with in homalg 
using the respective package. Further, and without any extra help from the ring package, 
homalg can automatically compute over residue class rings of any supported ring. 

• PIR [Bar07] is one more tiny package, or rather a pseudo-package, that makes 
Maple's built-in facilities for dealing with integers and some other principal ideal 
rings available to homalg. (Prime sub fields Q and Z/pZ and their finite field ex- 
tensions, realized as primitive extensions, rational function fields over the previous 
fields, the integers Z, the GAUSSian integers Z[\/— 1] and univariate polynomial 
rings Z/pZ[x], where p is a prime, Q[x] and K[x], where K is a rational function 
field over a finite extension of Q, realized as a primitive extension.) 

• Involutive [BCG + 03] implements the involutive basis technique of V. P. Gerdt 
and Y. A. Blinkov in Maple. An involutive basis is a special kind of Grobner 
basis for an ideal of a polynomial ring or, more generally, for a submodule of a free 
module over a polynomial ring. Involutive bases have nice combinatorial properties 
[PR05], and the algorithms designed by V. P. Gerdt and Y. A. Blinkov [Ger05, 
GB98a, GB98b] compute them efficiently. In fact, these algorithms provide an 
efficient alternative to Buchberger's algorithm [Buc06] to compute Grobner 
bases. Involutive restricts to particular involutive bases, namely Janet bases. It 
also provides an interface to a C++ implementation of the involutive basis technique 
which can be used to call the fast routines when needed as well as to switch to 
these fast routines for the whole Maple session. (Commutative polynomial rings: 
S[xi, . . . , x n ], where S is either Z or a field existing in Maple.) 

• Janet [BCG + 03] implements the involutive basis technique for computing Janet 
bases of linear systems of partial differential equations. (Differential algebras over 
differential fields: K[-t^-, . . . , gf-], where K is a differential field which exists in 
Maple.) 

• JanetOre [Rob06, Rob07] generalizes Involutive from commutative polynomial 
rings to certain iterated skew polynomial rings. In particular, it computes Janet 
bases for left ideals in Ore algebras [CS98]. (K[d;a, 5], where K is a polynomial 
ring over a field, d a new indeterminate, a is a certain automorphism of K and 5 a 
a-derivation of K, and iterated extensions of this kind.) 

• OreModules [CQR07] is a Maple package for the study of structural properties 
of linear systems over Ore algebras, i.e. linear equations involving certain linear 
functional operators which can be considered as elements of an Ore algebra. By 
default, it uses the Maple package Ore_algebra [CS98] to compute Grobner bases, 
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but these calls can also be switched to JanetOre. (Ore algebras [CS98] and the 
iterated skew polynomial rings from the previous point.) 

homalg is also able to make use of various normal form algorithms for modules resp. 
special types of modules over various rings, which are used to provide a standard form for 
a presentation of these modules: 

• PIR uses the Smith normal form for (Maple-built-in) principal ideal rings. 

• Janet optionally uses the JACOBSON normal form for univariate differential rings, 
i.e. rings of the form K[d], where K is a differential field with d a derivation of K. 

• Involutive optionally uses the extension package QuillenSuslin written by ANNA 
Fabianska [Fab07, FQ07] implementing the Quillen-Suslin theorem to compute 
a free basis of a projective module over a polynomial ring (which is then free by 
the theorem). A similar extension package is planned for OreModules. 

• OreModules optionally uses the extension package Stafford [QR] which computes 
a free basis for a stably free module of rank at least 2 over the Weyl algebras 
k[x±, . . . , x n , d±, . . . , d n ] and k(xi, . . . , x n )[dx, . . . , d n ], with k a field of characteristic 
0. 

Appendix B. The Maple implementation 

B.l. Presentations. A presentation of a module in the current Maple implementation of 
homalg is a list 19 containing as first entry the list of generators and as second entry the 
list of relations. The third entry is a string delimiter to optically indicate the end of the 
presentation. This string, unless changed by the user, defaults to "Presentation". The 
remaining entries provide extra information about the presented module, e.g. its Hilbert 
series in case the ring is the polynomial ring. This extra information can only be provided 
by the ring-specific package. 

For M = coker(M) = R lxl ° / R lxll Vl the list of the concrete generators are numbered by 
abstract generators being the l standard basis row vectors of the underlying free module 
R lxl °. The list of relations simply contains the rows of the matrix M. An illustration is 
given in Figure 1 of Example C.l. 

B.2. The morphism part of functors. The name convention for functors used in the 
Maple implementation of homalg is as follows: If the procedure implementing the object 
part of a functor is called F, then the procedure implementing the morphism part is called 
FMap. It is defined using the procedure FunctorMap applied to the object part procedure 
F. The several pieces of code for the morphism part of the functors Cokernel, Kernel, 
Def ectOfHoms and Hom_R, and the two bifunctors TensorProduct and Horn, reproduced on 
the web [BR07a], demonstrate how FunctorMap unifies the definition of the morphism part 
of all functors in homalg. 

In future implementations of homalg the two procedures implementing the object and 
morphism part of a functor will be unified in one. The unified procedure will be able to 



19 In future implementations we will make use of object oriented data structures, which encapsulate all 
this information. 
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recognize if it has been applied to an object, to a morphism or even to complexes. This 
will be an easy task, once we strictly use structures throughout homalg. 

B.3. Encapsulating functors. A functor is fully defined when both parts are defined, 
i.e. its object part and its morphism part. If the functor is a multi-functor, then sev- 
eral morphism parts have to be defined, homalg accesses all these parts of a functor via 
a so called encapsulation. It incorporates both parts and also takes care of the possi- 
ble multi-functoriality. The implementation of the encapsulations of the functors Hom_R, 
TensorProduct and Horn can be viewed under [BR07a]. 

B.4. Composition of functors. As was mentioned in Section 5, composing two functors 
is an easy task since one simply has to compose their actions on objects and on morphisms. 
The procedure responsible for this is called ComposeFunctors. The implementation of the 
composition of the functors HomHom_R and HomHom can be viewed under [BR07a]. 

B.5. Applying functors to complexes. As mentioned in Section 7, functors should be 
applied to complexes of modules, rather than to single modules. Thereby homalg again 
makes use of the encapsulation of functors. In this implementation, there are two pro- 
cedures depending on whether one is dealing with a covariant or a contravariant functor. 
They are called FunctorOnSeqs and Cof unctorOnSeqs. The implementation of the func- 
tors Horn and HomHom on complexes can be viewed under [BR07a]. 

B.6. Derivation of functors. The left derivation procedure for covariant functors is 
called Lef tDerivedFunctor. The faster derivation procedure for right exact functors, 
for example — <E> L, referred to at the end of Subsection 7.1 is called LeftDerived- 
RightExactFunctor. 

The right derivation procedure for contravariant functors is called RightDerivedCo- 
f unctor. The faster derivation procedure for left exact contravariant functors, for example 
Hom(— , L), referred to at the end of Subsection 7.1 is called RightDerivedLef tExact- 
Cof unctor. 

The implementation of the left derived functor LHomHom (of HomHom with respect to its 
first argument) and the right derived functor Ext (of Horn with respect to its first argument) 
can be viewed under [BR07a]. 

B.7. How FunctorMap works. Let F denote the object part procedure of a functor F. 
First FunctorMap asks F if the underlying functor is co- or contravariant. If F is homalg- 
basic it has to know the answer itself. If F is defined as a composition Fi o F 2 then the 
question is passed to the procedure ComposeFunctors, which decides the answer by asking 
F\ and F 2 (this is recursive) . If F is defined as a derivation L q G or R 9 G, using one of the 
derivation procedures, then F passes the question to the latter, which decides the answer 
by asking G (this is recursive) . Any recursion ends when a homalg-basic functor is reached. 

FunctorMap then asks F if it is defined using the procedure ComposeFunctors. If F is 
homalg-basic it ignores the question. If F is defined using one of the derivation procedures, 
the question is passed to the latter, which ignores it. If F is indeed defined as a composition 
FxoF 2 , then F passes the question to ComposeFunctors, which returns the two functors Fx 
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and F 2 in both their parts. FunctorMap can now easily construct the morphism part of F 
by composing the morphism parts of Fi and F 2 . In case F is not defined as a composition, 
FunctorMap asks it if is defined by derivation. If F is homalg-basic it ignores the question. 
If F is indeed defined as a derivation L q G or R q G, using one of the derivation procedures, 
then F passes the question to the latter, which returns the functor G in both its parts and 
a procedure based on ResolutionOf Seq to compute <p q out of if (cf. (Lift), p. 18). With 
these two ingredients FunctorMap is able to construct the morphism part of F as described 
in Subsection 7.2. If F is defined neither by composition nor by derivation, i.e. its is homalg- 
basic, then it is asked by FunctorMap to return 20 its hull functor Hull F together with the 
natural embedding (cf. (Hull), p. 10). Again this suffices to construct the morphism part 
of F as described in Section 5. 

Appendix C. Examples 

For further examples we refer to [BR06a, BR06b] and the site of homalg [BR07c]. 

A nice application is the tiny homalg-based package conley that computes C-connection 
matrices of graded module octahedra/braids of Morse decompositions in dynamical system 
theory [BR]. 

For the sake of demonstration we wrote a tiny package called alexander [BR07b], which 
relies on homalg and computes simplicial homology and cohomology. Future implementa- 
tions of homalg are planned to enable more serious applications to topology. 

C.l. Example 1. In this example we compute a module of homomorphisms. 

> restart : 

> with(Involutive) : with(homalg) : 

Specify the homalg-table of the ring package Involutive: 

> RPI :=' Involutive/homalg' ; 

RPI := Involutive / homalg 
Use the ring package Involutive as the default ring package: 

> ' homalg/def ault ' : =RPI ; 

homalg / default := Involutive/ homalg 
Define the ring R = Q[x, y, z\. 

> var : = [x , y , z] ; 

var := [x, y, z] 

> K:=Cokernel( [[x,y,0] , [x~2,y~2,0] , [x~3,y~3,z~3] ] ,var) ; 

K := [[[1, 0, 0] = [1, 0, 0], [0, 1, 0] = [0, 1, 0], [0, 0, 1] = [0, 0, 1]], 
[[x, y, 0], [0,xy- y 2 , 0], [0, 0, z% "Presentation", 

14 6 

3 + 8 s + 14 s 2 + s 3 ( + ) [Mj 6j Q]] 

1 — s 1 — sr 



'The technical details to realize this heavily depend on the implementation. 
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Figure 1 . A module of homomorphisms between two modules over R 
Q[x,y,z] with Involutive 



[[1,0,0] 



"10 0" 




" 




,[0,1,0] = 




1 






,[0,0, 1] 





x - y 



[[0,0,y],[x 2 - xy,0,-x],[0,z s ,0]}, 
"Presentation" , 
generators 3 + 8 s + 14 s 2 + s 3 ' 



+ 



(i-s) -T- jizrjyz 1 1 

relations [14, 6, 0] ] 

Hilbert series 

Cartan characters 



> L:=Cokernel( [[x,y]] ,var) ; 

L : = 
2 + s 



1, 0] = [1, 0], [0, 1] = [0, 1]], [[x, 
2 2 1 



"Presentation" . 



+ 



, [2, 2, 1]] 



•1-3 {1-8) 

Compute the module of homomorphisms Honifj(L, K) (see Figure 1): 
> hom : =Hom (L , K , var ) ; 



horn : = [[[1, 0, 0] 



3 + 8 s + 14 s 2 + s a 



[0, 1, 0] = 

"Presentation' 



1 
1 

[[0, 0, y], [x 2 -xy, 0, -x], [0, z , UJJ , 



-y 
x 



[0, 0, 1] 




-y + x 



C.2. Example 2. Let R := Q[x, y,z]. In this example we want to study two non- 
equivalent extensions 0— > K — > M — > L — > and ^ K — > N — > L with K a 
torsion module and L a torsion free module. Our goal is to use the notion of functor to 
reveal that these extensions are not only non-equivalent but also non-isomorphic. For this 
we define two functors Fm and Fn and study their behavior when applied to complexes. 
Concretely, we apply Fm resp. F^ to — * K — > M ^ L — > which we refer to as 

S : ^ K ^ M ^ L ^ 0. 

In the following we consider 

F Q = Ext(l, -,Q) = Ext^-, Q) = R 1 Hom fl (- Q) 



> restart ; 

> with(Involutive) : with(homalg) : 

Use the ring package Involutive as the default ring package: 
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> RPI :=' Involutive/homalg' ; 

RPI := Involutive / homalg 

> ' homalg/def ault ' : =RPI ; 

homalg / default := Involutive /homalg 
We force the Maple ring package Involutive to use the external program JB, which 
a C++ implementation of the involutive division algorithm: 

> InvolutiveOptions("C++") ; 

Define the ring R = Q[x, y, z\: 

> var : = [x , y , z] ; 

var := [x, y, z] 
The two presentation matrices M and N: 



> MM:=matrix( [[x*z, z*y, z~2, 0, 0, y] , [0, 0, 0, z~2*y-z~2, z~3, x*z] , 

> [0, 0, 0, z*y~2-z*y, z~2*y, x*y] , [0, 0, 0, x*z*y-x*z, x*z"2, x~2] , 

> [-x*y, -y~2, -z*y, x~2*y-x~2-y+l , x~2*z-z, 0], [x~2*y-x~2, x*y~2-x*y, 

> x*z*y-x*z, -y~3+2*y~2-y, -z*y~2+z*y, 0]]); 



MM :-- 



X z 


zy 


z 2 










11 











z 2 y — z 2 


z 3 




X z 











zy 2 - zy 


z 2 y 




xy 











xzy — x z 


2 

X Z 




x 2 


-xy 


-y 2 


-zy 


x 2 y — x 2 — y + 1 


2 

X z — 


z 





2 

y — x 


xy 2 — xy 


x zy — x z 


-y 3 + 2y 2 -y 


-zy 2 + 


zy 






> 
> 
> 
> 



NN:=matrix( [[x*y, y~2, z*y, 0, 0, z] , [0, 0, 0, z*y~2-z*y, z~2*y, 
x*z] , [x~2*z, x*z*y, x*z"2, -z"2*y+z"2, -z~3, 0], [0, 0, 0, 
y~3-2*y"2+y, z*y"2-z*y, x*y-x] , [0, 0, 0, x~2*y-x~2-y+l , x~2*z-z, y] , 
[x~3, x"2*y, x"2*z, -x*z*y+x*z, -x*z~2, 0]]); 



NN :-- 



xy 


ir 


zy 








Z 











zy 2 - zy 


z 2 y 


X z 


2 

X z 


xzy 


2 

X Z 


—z 2 y + z 2 


-z 3 














y 3 - 2 y 2 + y 


zy 2 - zy 


xy - 











x 2 y — x 2 — y + 1 


x 2 z — z 


V 


X 3 


2 

x y 


2 

X z 


—x zy + x z 


2 

—x Z 






X 



The two modules M := coker ( R 



>lx6 



R 



1x6 



and iV := coker f R 



>lx6 



R 



1x6 



> M : =Cokernel (MM , var) ; 
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M := [[[1, 0, 0, 0, 0, 0] = [1, 0, 0, 0, 0, 0], [0, 1, 0, 0, 0, 0] = [0, 1, 0, 0, 0, 0], 

[0, 0, 1, 0, 0, 0] = [0, 0, 1, 0, 0, 0], [0, 0, 0, 1, 0, 0] = [0, 0, 0, 1, 0, 0], 

[0, 0, 0, 0, 1, 0] = [0, 0, 0, 0, 1, 0], [0, 0, 0, 0, 0, 1] = [0, 0, 0, 0, 0, 1]], [ 

[xz, zy, z 2 , 0, 0, y], [0, 0, 0, z 2 y - z 2 , z 3 , xz], [0, 0, 0, zy 2 - zy, z 2 y, xy], 

[0, 0, 0, x zy — xz, x z 2 , x 2 ], [x 2 z, x zy, xz 2 , 0, 0, xy], 

[-x y, -y 2 , -z y, x 2 y - x 2 - y + 1, x 2 z - z, 0], 

[x 2 y — x 2 , xy 2 — xy, x zy — xz, —y 3 + 2y 2 — y, —zy 2 + zy, 0], 

[0, 0, 0, z y — z, z 2 , x 3 — y 2 }}, "Presentation", 

3 s 2 2 s 2 3 2 s 2 „ 1 s 

1 1 1 1 1- s 4- s -I 1 

1-5 1-5 (1-S) 2 (1-S) 3 (1-S) 2 1-5 (1-s) 2 ' 

[34, 14, 3]] 

> N : =Cokernel (NN , var ) ; 

N := [[[1, 0, 0, 0, 0, 0] = [1, 0, 0, 0, 0, 0], [0, 1, 0, 0, 0, 0] = [0, 1, 0, 0, 0, 0], 

[0, 0, 1, 0, 0, 0] = [0, 0, 1, 0, 0, 0], [0, 0, 0, 1, 0, 0] = [0, 0, 0, 1, 0, 0], 

[0, 0, 0, 0, 1, 0] = [0, 0, 0, 0, 1, 0], [0, 0, 0, 0, 0, 1] = [0, 0, 0, 0, 0, 1]], [ 

[xy, y 2 , zy, 0, 0, z], [0, 0, 0, zy 2 - zy, z 2 y, x z], [x 2 z, xzy, xz 2 , -z 2 y + z 2 , -z 3 , 

[0, 0, 0, y 3 -2y 2 + y, zy 2 - zy, xy - x], [0, 0, 0, x 2 y - x 2 - y + 1, x 2 z - z, y], 

[x 2 y, xy 2 , xzy, 0, 0, x z], [x 3 , x 2 y, x 2 z, —xzy + xz, —xz 2 , 0], 

[0, 0, 0, xy 2 z — xzy, xy z 2 , x 2 z], [0, 0, 0, 0, 0, x 3 z — zy 2 — x z], 

[0, 0, 0, xy 3 — 2 xy 2 + xy, xy 2 z — xzy, x 2 y — x 2 ], 

3 s 2 

[0, 0, 0, 0, 0, x 3 y - x 3 - y 3 - xy + y 2 + x]\, "Presentation", h -, + 2 s 2 

1 — s (1 — s) 2 

3 2 s 2 s 2 „ 1 s 3 s 

The torsion submodule K := t(M) of M (and N): 

> K : =TorsionSubmodule (M , var) ; 

K := [[[1, 0, 0] = [0, 0, 0, 0, 0, 1], [0, 1, 0] = [0, 0, 0, y - 1, z, 0], [0, 0, 1] = [x, y, z, 0, 0, 
], [[y, 0, z), [0, -zy, xz], [xz, z 2 , 0], [0, -y 2 + y, xy - x], [xy, zy, 0], [0, x 2 - 1, -y), 

Is 2 s 

[x 2 , xz, 0]], "Presentation", s + + + + , [7, 3, 0]] 

1 — s (1 — s) 2 (1 — s) 2 1 — s 

The embedding — > K — 1 -+ M: 



> alphal :=TorsionSubmoduleEmb(M,var) ; 
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1 
al:= OOOy-lzO 
x y z 
The torsion free factor L := M/K = M/t(M) of M (and N): 

> L:=Cokernel(alphal,M,var) ; 

L := [[[1, 0, 0, 0, 0] = [1, 0, 0, 0, 0, 0], [0, 1, 0, 0, 0] = [0, 1, 0, 0, 0, 0], 
[0, 0, 1, 0, 0] = [0, 0, 1, 0, 0, 0], [0, 0, 0, 1, 0] = [0, 0, 0, 1, 0, 0], 
[0, 0, 0, 0, 1] = [0, 0, 0, 0, 1, 0]], [[0, 0, 0, y - 1, z], [x, y, z, 0, 0]], "Presentation' 1 
2 3 

The natural epimorphism M — ^ L — > 0: 

> alpha2 : =CokernelEpi (alphal ,M,var) ; 



a2 :-- 



The short sequence S : — > K — l -> 

> IsShortExactSeq(K, alphal, M,alpha2,L, var) ; 

true 

The functor Ext}j(— , N) applied to the short exact sequence S: 

> coseqN : =Ext0nSeqs ( 1 , N , [alpha2 , alphal] , [L , M , K] , var ) : 

This says that the map Ext^(L, N) — 2 ' ; ExtJj(M, N) is not injective, i.e. that the 

8° i 

connecting homomorphism H.omn(K, N) — > Ext fi (L, N) is non-trivial: 

> IsShortExactSeq(op(MakeCoseq(coseqN)) , var , "VERBOSE") ; 





' 1 








o o ■ 









1 


















1 


















1 















1 






. 








o o . 






L 




is exact 



'horns" = true, "cmps" = true, "defs" 


1 



z, y, x], "Presentation", 1, [0, 0, 0]], true, true] 



The functor Ext^(— , M) applied to the short exact sequence S: 

> coseqM : =Ext0nSeqs ( 1 , M , [alpha2 , alphal] , [L , M , K] , var) : 

Whereas the resulting sequence Ext^S 1 , M) is again exact. In particular, the connecting 

8° i 

homomorphism Homji(K, M) — > Ext i? (L, M) is trivial, i.e. the zero map: 

> IsShortExactSeq(op(MakeCoseq(coseqM)) , var , "VERBOSE") ; 

true 
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C.3. Example 3. This example demonstrates by a computation over a noncommutative 
Ore domain how one can use homalg to solve a linear system of differential-difference 
equations: 

> restart : 

> with( JanetOre) : with(homalg) : 

Specify the homalg-table of the ring package JanetOre: 

> RP0:=' JanetOre/homalg' ; 

RPO := JanetOre/homalg 
Use the ring package JanetOre as the default ring package: 

> ' homalg/def ault ' : =RP0 ; 

homalg / default := JanetOre/homalg 
Ask homalg to check if a stably free module with a finite free resolution 21 is free, cf. [QR, 
Remark 51]: 

> homalg_options("check_rank_l"=true) ; 

"JanetOre/homalg", "check_if_stably_free_rank_l_FFRJs_free" = true 
Define the Ore domain R = Q[t) [D, t ^ t, |] [5, t \-> t + 1, 0]: 

> Ore: = [[t,D,delta] , [] , [weyl(D.t) .shift (delta, t)]] ; 

Ore := [[t, D, 5], [], [weyl(D, t), shift t)}} 
Find the general solution to the following linear system of differential-difference equa- 
tions: 

u{t) + u{t + 2) + D(v)(t + l) -v(t + 2) + w{t + l) = 0, 
tD (u) (t + 1) + u (t + 3) - 2 tu (t + 1) + 2 D (u) (t + 1) - 2 u (t + 1) 

-v{t + 3) + 2v(t + 2) + w(t + 2) = 0. 

This system can be written as the following matrix operators A applied to the section 
( u(t) v(t) w(t) ) . (Caution: A monomial in t, D and 5 must be read as a monomial 
with powers of t on the left of a monomial in the commuting D and 5: read in the following 
both tD, Dt as tD, and both t5, St as tS): 

> A := matrix( [[l+delta~2, delta*D-delta~2 , delta], 

> [t*delta*D+delta~3-2*t*delta+2*delta*D-2*delta, -delta"3+2*delta~2, 

> delta~2]]); 

[ 1 + 5 2 5D- 5 2 5 ' 

[ t5D + 5 3 -2t5+25D -25 -5 3 + 25 2 5 2 
Studying the system is equivalent to studying the cokernel M of A e R mxn viewed as a 
morphism R — > R . This is based on the following observation: Let 5F be a function 
space, where one wants to search for solutions of the system, then Hom^(coker(A), £F) = 
Solj-(yl) := {rj G 3 rnxl | A?7 = 0}, where n is the number of columns of A. 



Without a finite free resolution homalg cannot check stably freeness, cf. [QR, Remark 28]. 



homalg 
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M := [[[1, 0, 0] = [1, 0, 0], [0, 1, 0] = [0, 1, 0], [0, 0, 1] = [0, 0, 1]], 
[[1 + 5 2 , 5B-5 2 , 5], [t5B + 5 3 -2t5 + 25B-25, -5 3 + 25 2 , 5 2 ]}, 

17 8 1 . . 
+ + [17, 8, 1]] 



"Presentation" ,3 + 9s + 17s 2 + s s 



•1-5 (1-S) 2 (1-S) 

The torsion free factor F of M turns out to be free of rank 1 (the above mentioned 
option "check_rank_l" was used): 

> F : =TorsionFreeFactor (M , Ore) ; 

1 



F := [[1 = [-5, 5 - D, -1]], [0], "Presentation" 

The natural epimorphsim M —> F: 
> nu:=TorsionFreeFactorEpi(M,Ore) ; 



1-3 



[0, 0, 1]] 



v :- 



S 
t 



-2 + 5 + t5-tD-5 2 j 
Since F is free one can use Leftinverse to compute a split M F (cf. 3.1.3,(1)): 

> chi :=Leftinverse(M,nu,F,0re) ; 

X :=[-S 6-T> -1 ] 
Now compute the torsion submodule T of M: 

> T : =TorsionSubmodule (M , Ore) ; 

T := [[[1, 0] = [t+1, -5, 0], [0, 1] = [-2 5, -tS + tD + 5 2 + l, t}}, 
[[1 + 5 2 , 8], [-2 + D, -25 + 5 D], [-2 5 + SB, 0]], "Presentation", 
6 3 



2 + 6 s + s 2 



[6, 3, 0]] 



1-s (l-s) 2 ' 
And the natural embedding T —> M: 

> iota: =TorsionSubmoduleEmb(M, Ore) ; 

t+1 -5 

-25 -t5 + tB + 5 2 + 1 t 
After a few tests we found an element fi G T 

> mu: = [ [delta, 1] ] ; 

A* := 1]] 

which turned out to be a cyclic generator of T (below one identifies [i with the morphism 

R lXl _^ T . fi ); 

> IsSurjective(mu,T,Ore) ; 

true 
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H is nothing but T rewritten on this cyclic generator. It further turns out that the cyclic 
generator \x of the torsion submodule T satisfies a single simple relation S 2 (D — 2)fi = 0: 

> H : =Image (mu , T , Ore) ; 

H := [[1 = [t5, tB-t5 + l, t\], [-2 5 2 + S 2 D], "Presentation", 

1 + 3 s + 6 s 2 + s 3 (-^- + 3 - ), [6, 3, 0]] 
1 — s (I — s) z 

e is nothing but t rewritten on the cyclic generator fi: 

> epsilon:=ComposeMaps(mu, iota, M, Ore) ; 

e:=[tS tD-tS+1 t] 

> IsHom(H,epsilon,M,Ore) ; 

true 

> IsInjective(H,epsilon,M,Ore) ; 

true 

Now construct the isomorphism a from the direct sum S := F © H onto M: 

> alpha :=RP0 [matrix] (RPO [UnionOf Rows] (chi , epsilon) ) ; 

_ [ -5 5-D -1 " 
a: ~ [t5 tD-tS+1 t 

> S:=DirectSum(F,H,Ore) ; 

S := [[[1, 0] = [1, 0], [0, 1] = [0, 1]], [[0, -2 5 2 + 5 2 B\], "Presentation" , 

2 + 6 s + 12 s 2 + + + (T^ja)' t 12 ' 6 ' !]] 

> IsHom(S, alpha, M, Ore) ; 

true 

> IsBijective(S, alpha, M, Ore) ; 

true 

N is now nothing but M rewritten on the images of the two generators of S. So M (or 
N) is a module generated by two generators, one free and one subject to a single simple 
relation: 

> N :=Image (alpha, M, Ore) ; 

N:= [[[1, 0] = [-6, 5-D, -1], [0, 1] = [t6, tD-tS + 1, t}}, [[0, —2 5 2 + 5 2 D]], 

12 6 1 

"Presentation", 2 + 6 s + 12 s 2 + s 3 (- + - + [12, 6, 1]] 

1 — s (1 — s) z (1 — s)- 1 

Using sections, the above presentation reads: The original system on the three unkown 

functions u(t), v(t) and w(t) is equivalent to a system on two unknown functions f(t) and 

g(t), where only g(t) satisfies the single simple equation D(g)(t + 1) — g(t + 1) = 0. It is 

also given how to express / and g in terms of u, v and w. 
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> PresentationOnSections(N,Ore [1] ,0re [3] , [u,v,w] , [f ,g] ) ; 

[[f(t) = -u(t + 1) + v(f + 1) - D(v)(t) - w(t), 

g(t) =tu{t + l)+tD(v){t)-tv(t + l)+v(t) + tw(t)\, [-2g(t + 2) + D(#)(t + 2)], 

12 6 1 , . 
+ + [12, 6, 1]] 
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1-s (1-s) 2 (1-s) 
/(£) is therefore a free function and — C exp(t), where C is an arbitrary constant: 

> sol := [f =unapply(f (t) ,t) , g=unapply(C*exp(2*t) ,t)] ; 

sol:=[f = {t^l{t)),g={t^Ce^)] 
In order to solve the original system one needs to express u, v and w in terms of / and 
g. To this end use Leftinverse to compute r\ := a -1 (cf. 3.1.3,(2)): 

> eta :=Leftinverse(S, alpha, M, Ore) ; 

5 
r) := t 1 

-2 + 6 + t5-tD- 5 2 5 -T> 
Again, L is nothing but S (= M) rewritten on the images (under 77) of the three original 
generators of M: 

> L : =Image (eta , S , Ore , "USE_IMAGE_OF_GENERATORS" ) ; 

L ■= [[[1, 0, 0] = [5, 0], [0, 1, 0] = [t, 1], [0, 0, 1] = [-2 + d + td-tD-d 2 , 5-D}}, 
[[1+S 2 , 5D-5 2 , 8], [t5D + 5 3 -2t5 + 25D- 25, -5 3 + 25 2 , 5% 

"Presentation", 3 + 9s + 17s 2 + s 3 (-^— + - 8 N9 + - 1 - ), [17, 8, 1]] 

1 — s (1 — s) z (1 — s) 6 

And again, using the JanetOre procedure PresentationOnSections, the above presen- 
tation is rewritten: 

> P :=PresentationOnSections(L , Ore [1] , Ore [3] , [f ,g] , [u,v, w] ) ; 

P:=[[u(t)=f(t + l),v(t)=tf(t)+g(t), 

w(t) = -2 f(t) +{(t + l)+ t f(t + 1) - t D(/)(t) - f(t + 2) + g(t + 1) - D(g)(t)}, [ 
u(t) + u(t + 2) + D(v)(t + 1) - v(t + 2) + w(t + 1), t D(u)(t + 1) + u(t + 3) 
- 2 1 u(t + 1) + 2 D(n) (t + 1) - 2 u(t + 1) - v(t + 3) + 2 v(t + 2) + w(t + 2)], 

17 8 1 



"Presentation" , 3 + 9 s + 17 s 2 + s £ 



:i-s) 



r)» [17, 8, 1]] 



1-s (1-s) 2 
The most general C 00 (IR)-solution is given by: 

> Sol : =eval (RP0 [matrix] (GeneratorsOf Presentation(P) ) ,sol) ; 

f(t + 1) 

Sol:= tf{t) + CeW 

-2 f(t) + f(t + 1) + 1 f(t + 1) - t D(/)(t) - f(t + 2) + C e( 2 *+ 2 ) - 2 C e^ 2 *) 
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Besides, one has the following uniqueness property: If / and C lead to the same solution 
then /= / and C = C. 

Test the general solution with the JanetOre procedure JApplyMatrix: 

> JApplyMatrix (A, Sol, Ore [1] ,0re[3]) ; 

" " 
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